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1.0  INTRODUCTION 


The  effect  of  unsteady  motion  of  an  airfoil  on  Its  stall  behavior  Is  of 
considerable  Interest  to  many  practical  applications  Including  the  blades  of 
helicopter  rotors  and  of  axial  compressors  and  turbines^.  Experiments  with 
oscillating  airfoils,  for  example,  have  shown  that  the  flow  can  remain  attached 
for  angles  of  attack  greater  than  those  which  would  cause  stall  to  occur  In  a 
stationary  system.  This  result  appears  to  stem  from  the  formation  of  a  vortex 
close  to  the  surface  of  the  airfoil  which  continues  to  provide  lift.  It  Is 
also  evident  that  the  onset  of  dynamic  stall  depends  strongly  on  the  airfoil 
section  and  as  a  result  great  care  Is  required  In  the  development  of  a 
calculation  method  which  will  accurately  predict  this  behavior. 

In  principle,  the  prediction  of  dynamic  stall  can  be  accomplished  by  solving 
the  Reynolds-averaged  Navler-Stokes  equations  or  their  reduced  forms.  A 
turbulence  model  Is  required  and  Is  presumed,  with  reasonable  supporting 
evidence,  see  Ref.  2,  to  be  uninfluenced  by  the  Imposed  unsteadiness.  Several 
papers  have  been  prepared  with  calculations  of  this  type,  see  for  example  Ref. 
3,  and  Involve  the  solution  of  equations  with  two  diffusion  terms  as  well  as 
parabolized  forms  and  thin-layer  approximations.  An  alternative  Is  to  make 
use  of  Interactive  boundary-layer  theory  whereby  Invlscld  and  boundary-layer 
equations  are  solved  and  allowed  to  Influence  each  other  by  an  Iterative 
scheme. 

Extensive  Investigations  with  an  Interactive  approach  have  been  reported  by 

4 

Cebecl  et  al.  and  show  that  the  Incompressible  flow  and  performance 
characteristics  of  airfoils  can  be  predicted  accurately  and  efficiently  for 
high  and  low  Reynolds  numbers  and  for  a  range  of  angles  of  attack  up  to  and 
Including  stall.  At  Incidence  angles  higher  than  stall,  however,  this 
procedure,  was  unable  to  predict  the  airfoil  performance  due  to  relatively 
large  regions  of  flow  separation  on  the  surface  and  In  the  wake.  Near  stall, 
the  value  of  the  trailing  edge  displacement  thickness  approached  10%  of  the 
chord  and  the  numerical  method  could  not  provide  converged  solutions.  The 
predictions  of  this  Interactive  boundary-layer  are  similar  to  those  obtained 
from  solutions  of  thin-layer  Navler-Stokes  by  the  ARC-20  method^’6  for 
angles  of  attack  up  to  and  Including  stall.  Fig.  1.  It  has  been  shown  In  Ref. 

4  that  the  Interactive  flow  calculations  without  the  wake  effect  and  for 
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angles  of  attack  greater  than  that  of  stall,  yielded  lift  coefficients  which 
Increased  with  Incidence  angle  almost  In  the  same  way  as  those  computed  with 
the  thin-layer  Navler-Stokes  approach  with  the  wake  effect  Included.  When  the 
wake  effect  was  Included  In  the  Interactive  boundary- layer  calculations ,  the 
results  agreed  more  closely  with  measurements  but  could  not  be  extended  beyond 
the  stall  angle. 

More  recently,  the  Interactive  method  has  been  Improved  to  permit  calculations 
for  angles  of  attack  greater  than  that  of  stall  and  the  results  have  been 
shown  to  have  the  correct  behavior.  To  achieve  this,  modifications  were  made 
to  the  Iterative  procedure  and  to  the  method  of  calculating  the  wake.  These 
Improvements  are  described  In  Ref.  8  and  are  necessary  where  results  are 
required  at  angles  of  attack  corresponding  to  stall  and  post  stall. 

The  Interactive  boundary-layer  method  for  steady  flows  has  also  been  extended 

to  study  the  laminar  separation  and  reattachment  near  the  leading  edge  of  a 

9 

thin  oscillating  airfoil  ,  but  the  calculation  of  flow  over  practical 
airfoils  Involves  laminar  and  turbulent  flows  and  the  Inclusion  of  the 
upstream  Influence  of  the  wake  requires  careful  step-by-step  development  and 
evaluation,  as  has  been  done  for  steady  flows.  The  use  of  a  quasi-steady 
approach  to  unsteady  flows  thus  represents  an  essential  building  block  In  a 
progression  towards  an  Interactive  calculation  method  which  solves  unsteady 
equations  even  though  the  latter  Is  likely  to  represent  a  much  wider  range  of 
oscillation  frequencies.  The  extent  of  the  differences  can  be  quantified  only 
by  comparing  results  from  both.  The  present  study  permits  assessment  of  a 
quasi-steady  method  over  a  range  of  angles  of  attack  and  frequency  In  terms  of 
convenience,  accuracy,  and  the  computational  cost.  It  should  be  remembered 
that  the  main  purpose  of  boundary-layer  calculations  Is  to  Include  the  viscous 
effects  In  the  Invlscld  method  and,  since  this  Is  achieved  In  an  Iterative 
manner.  It  Is  useful  to  use  a  quasi-steady  approach  In  the  early  stages  of  the 
unsteady  flow  calculations.  This  strategy  can  Improve  the  efficiency  of  the 
calculations,  reduce  the  computer  time  and  avoid  possible  difficulties  that 
may  arise  when  the  unsteady  boundary-layer  equations  are  solved. 

The  calculation  of  the  wake  flow  becomes  more  difficult  with  Increasing  angle 
of  attack  and  Investigations  of  steady  flows  Indicate  that  near  and  past 
stall,  the  behavior  of  the  velocity  profile  at  the  trailing  edge  resembles 
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that  of  a  mixing  layer  with  a  huge  backflow.  The  thickness  of  the  upper 

velocity  profile  Is  an  order  of  magnitude  larger  than  the  lower  velocity 

profile,  and  the  upper  velocity  profile  contains  flow  reversal  which  Increases 

with  angle  of  attack,  whereas  the  lower  surface  velocity  profile  Is  almost 

lamlnar-llke  and  has  no  flow  reversal.  The  FLARE  approximation  In  which  the 

convective  term,  u  3u/3x,  In  the  x-momentum  equation  Is  set  equal  to  zero 

In  regions  of  negative  velocity,  has  yet  to  be  evaluated  for  wakes  of  this 

type  and,  Indeed,  the  only  real  test  of  the  FLARE  approximation,  first 

10  7 

suggested  by  Reyhner  and  Flugge-Lotz  ,  Is  that  of  Cebecl  et  al.  for 
steady  flows.  For  these  reasons,  a  useful  additional  step  towards  the 
long-term  objective  of  fully  unsteady-flow  calculations  Is  to  solve  the 
boundary- layer  equations  In  time-dependent  form  on  the  airfoil,  and  In 
quasi-steady  form  In  the  wake. 

The  purpose  of  this  paper  Is  to  present  calculations  of  unsteady  flows, 
obtained  with  the  quasi-steady  approach  and  to  assess  the  range  of  conditions 
for  which  this  approach  can  be  expected  to  provide  useful  results.  The 
Interactive  boundary-layer  method  Is  described  In  the  following  section  and 
the  results  presented  and  discussed  In  Section  3.  The  paper  ends  with  a 
summary  of  the  more  Important  conclusions. 
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2.0  INTERACTIVE  BOUNDARY-LAYER  METHOD 

The  unsteady  Interactive  boundary-layer  method  makes  use  of  the  panel  method 

developed  by  Platzer  and  his  student,  Teng  and  the  Inverse 

4 

boundary- layer  method  of  Cebecl  et  al.  .  The  solution  of  the  Invlscld  flow 

12 

equations  for  the  airfoil  Is  similar  to  that  of  Hess  and  Smith  and  uses 

1 3 

the  procedure  of  Basu  and  Hancock  for  the  wake.  The  shape  of  the  airfoil 
Is  represented  by  a  large  number  of  straight-line  segments  called  panels  (an 
Inscribed  polygon),  each  of  which  has  a  different  constant  source  strength. 

The  wake  Is  represented  by  a  series  of  free  vortices  shed  from  the  trailing 
edge,  each  having  a  constant  vortlclty  strength  whose  variation  along  the  wake 
Is  determined  by  the  time  history  of  the  circulation  about  the  airfoil. 


The  boundary- layer  method  Is  based  on  the  solution  of  the  two-dimensional 
steady  boundary-layer  equations  expressed  In  terms  of  an  eddy-viscosity, 
c  ,  so  that  the  continuity  and  momentum  equations  have  the  form, 


3u 

3x  + 


3u 

3x 


+  v 


3u 

3y 


=  u 


3u 
_ e 

e  3x 


3_ 

3y 


t?1 


(1) 

(2) 


and  are  subject  to  boundary  conditions 


y  =  0,  u  =  v  =  0,  y  -*  ®,  u  -»  u  (x)  (3) 

e 

which  Include  the  flow  on  the  airfoil.  With  y  =  0  now  denoting  the  dividing 
streamline  that  separates  the  upper  and  lower  parts  of  the  Invlscld  flow  In 
the  wake,  they  can  be  written  as 


y  -»  ®,  u  -*  u  (x) ; 

e 


y  =  0,  v  =  0; 


u  (x) 
e 


(4) 


In  order  to  avoid  the  breakdown  of  the  solutions  at  flow  separation,  the  outer 
boundary  condition  expresses  ug(x)  as  the  sum  of  Invlscld  velocity  u°(x) 

and  a  perturbation  velocity  <5u  (x)  computed  from  the  Hilbert  Integral 

e 

r  #  i  1  ,  b  d  i  r  x  .  (Iff 

Jue<x)  ■;  I  s  <“e‘  > 

xa 
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with  the  Interaction  scheme  confined  to  (x  ,  x.).  Introduction  of  a  dls- 

a  D 

Crete  approximation  to  this  Integral  enables  the  perturbation  velocity  to  be 


expressed  In  terms  of  the  geometric  coefficients  of  the  airfoil,  c  ,  and 
enables  the  edge  boundary  condition  to  be  written  In  a  form  In  which  a 


relationship  Is  provided  with  the  external  velocity  and  displacement 
thickness,  that  Is, 


Ue(*>  *  U°(x) 


J?,  C1j(Ue4*)j 


and  can  be  generalized  to  the  form 


ue(x)  -  ueK(x)  ♦  ^  c  1  j ( D j 


(5a) 


(5b) 


where  D  =  ug6*  and  ug  corresponds  to  the  Invlscld  velocity  distribution  which 
contains  the  displacement  thickness  effect  (6*)  computed  from  a  previous 

sweep,  as  discussed  next.  This  expression  was  suggested  by  Veldman  for  steady 
14 

flows  ,  and  can  also  be  viewed  as  an  empirical  formula  that  provides  a 
better  approximation  to  the  link  between  Invlscld  and  viscous  flow  equations 
than  the  relaxation  formula  used  by  teBalleur15  and  Carter16. 


The  solution  of  the  system  given  by  Eqs . ( 1  - 4 )  with  the  outer  edge  specified  by 
Eq.  (5)  and  with  the  algebraic  eddy  viscosity  of  Cebecl  and  Smith  Is  obtained 
by  Keller's  two-point  finite-difference  method  described  In  Ref.  4.  For  a 
given  pressure  distribution  and  chord  Reynolds  number,  R  (E  u  c/v),  the 

C  00 

boundary-layer  calculations  start  at  the  stagnation  point  and  are  performed 
for  a  laminar  flow  in  the  direct  mode  for  a  few  x-statlons  after  which  they 
are  performed  In  the  Inverse  mode.  Turbulent  flow  calculations  begin  at  the 
transition  location,  which  Is  specified  at  the  pressure  peak  or  computed  from 
formulas  such  as  that  of  Michel,17  and  extend  Into  the  wake.  The  calculated 
displacement  thickness  and  external  velocity  distribution  are  used  to 
determine  a  blowing  velocity  distribution,  vn  =  dO/ds,  so  that  the  Invlscld 
flow  equations  are  solved  again  subject  to  a  new  boundary  condition.  The 

Invlscld  velocity  distribution  Is  evaluated  directly  on  the  displacement 
surface  while  the  blowing  velocity  Is  applied  on  the  original  airfoil 
surface.  The  Kutta  condition  Is  also  satisfied  at  the  displacement  surface 
and  this  process  Is  repeated  on  an  Iterative  basis  until  the  solutions  of  both 
boundary-layer  and  Invlscld-f low  equations  converge. 
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This  procedure  Is  appropriate  as  long  as  there  Is  no  separation  on  the 
airfoil.  With  separation,  the  FLARE  approximation  was  used  and  proved  to  be 
satisfactory  when  the  separation  region  was  small.  As  the  extent  of  the 
separation  region  on  the  airfoil  increased  with  Increasing  angle  of  attack,  an 

4 

additional  Iterative  scheme,  based  on  a  homotopy  continuation  method  was 
Introduced  at  the  start  of  the  wake  calculations.  Thus,  an  Initial  velocity 
prof  1 le  was  defined  by 


u 


ref 


n(  u 

a 


“ref* 


(6) 


and  the  boundary- layer  equations  solved  at  the  first  x-statlon  on  the  wake  with 

n  =  0.  Here  u  ,  corresponds  to  a  nonseparating  velocity  profile  constructed 
ref 

arbitrarily  from  the  separated  velocity  profile  at  the  previous  station  that 
has  the  velocity  profile  denoted  by  u  .  The  calculation  was  then  repeated 
for  different  values  until  the  solutions  converged.  This  procedure  was 
applied  for  each  velocity  profile  In  the  wake  with  separation  and  was 
necessary  when  the  extent  of  the  tral 1 Ing-edge  separation  was  around  10%  of 
the  chord;  It  resulted  In  successful  calculations  of  the  lift  and  drag 
coefficients  of  several  airfoils  up  to  and  Including  stall.  With  the 
Improvements  described  In  Ref.  8,  It  can  now  be  extended  to  compute  post-stall 
flows  on  airfoils. 
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3.0  RESULTS  AND  DISCUSSION 


The  method  described  In  the  previous  section.  In  which  the  Invlscld  flowfleld 
Is  computed  by  the  .steady  panel  method  of  Ref.  11  and  the  viscous  flowfleld 
by  the  steady  Inverse  boundary-layer  method  of  Ref.  4,  was  used  to  compute 
flow  over  a  Sikorsky  SSC-A09  airfoil  subject  to  a  harmonic  oscillation  and 
ramp-type  motion  with  constant  pitch  rate.  At  first,  calculations  were 
performed  for  steady  flows  so  that  the  results  could  be  compared  with 
experimental  data  for  a  range  of  angles  of  attack  Including  post-stall,  as 
described  In  Subsection  3.1.  Next,  calculations  described  In  Subsection  3.2 
were  performed  with  the  quasi-steady  approach  by  Interacting  the  unsteady-flow 
panel  method  with  the  Inverse  steady  boundary-layer  method.  This  step  removed 
the  potential  difficulties  of  the  unsteady- flow  boundary  layer  and  allowed 
emphasis  on  the  unsteady- f low  panel  method.  The  quasi-steady  model  should,  of 
course,  be  able  to  simulate  flows  over  airfoils  executing  slow  motion  and  will 
overlap  the  final  work  of  unsteady  flow  calculations.  The  conditions  for 
which  quasi-steady  flow  calculations  are  appropriate  are: 


3u 
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3.1  Steady  Flow 

Calculations  were  performed  for  steady  flow  over  the  Sikorsky  airfoil  at  a 

6 

chord  Reynolds  number  of  2  x  10  with  the  location  of  transition  computed 

from  Michel's  formula,  except  where  the  boundary-layer  separated  upstream  of 

this  location.  In  which  case  transition  was  assumed  to  occur  at  the  separation 

point.  At  higher  angles  of  attack,  separation  occurred  at  a  very  short 

distance  from  the  pressure  peak  and  In  those  cases,  transition  was  assumed  to 

occur  at  a  very  small  distance  upstream  of  the  pressure  peak.  This  procedure 

allowed  the  calculations  to  be  performed  for  a  range  of  angles  of  attack, 

Including  post-stall,  and  to  be  compared  with  the  measurements  of  Lorber  and 
„  .  18 
Carta. 

Figures  2,  3  and  4  show  the  variations  of  lift  coefficient,  pressure 
coefficient  and  boundary-layer  parameters  as  a  function  of  angle  of  attack  and 
confirm  that  the  Interactive  method  Is  able  to  represent  measurements  for 
angles  up  to  and  beyond  that  of  stall.  The  close  agreement  between  the 
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calculated  lift  coefficients  and  the  measurements  of  Ref.  18  are  In  contrast 
to  the  solutions  of  the  Invlsc id-flow  equations  which,  as  expected,  Increase 
linearly  with  angle  and  do  not  recognize  that  the  airfoil  has  stalled.  The 
pressure- coefficient  distributions  of  Fig.  3  reveal  the  same  features  with  the 
Invlscld  and  Interactive  calculations  providing  similar  results  at  angles  of 
attack  below  II  degrees  and  Increasingly  different  results  as  the  angle  of 
attack  Increases;  again,  the  Interactive  results  are  In  close  accord  with 
measurements  except  perhaps  at  17  degrees  where  the  large  region  of 
upper-surface  separation  Introduces  Inaccuracies  In  measurements  and 
calculations. 


The  calculated  boundary-layer  properties  of  Fig.  4  show  the  expected  large 

variations  as  angle  of  attack  Increases  to  and  beyond  stall.  The  variations 

of  u  /u  are  strongly  dependent  on  the  angle  of  attack,  particularly  In 
0  00 

the  tral ling-edge  region  and  correctly  tend  to  unity  with  Increasing  x.  The 
displacement  thickness  on  the  upper  surface  achieves  very  high  values  and  on 
the  airfoil  Is  always  very  much  larger  than  that  on  the  lower  surface.  There 
Is  also  a  tendency  for  the  location  of  the  displacement-thickness  maximum  to 
move  from  the  trailing  edge  downstream  with  angle  of  attack  and  for  17  degrees 
It  occurs  at  x  =  1.07.  The  skln-f rlction  coefficient  confirms  that  the 
location  of  boundary- layer  separation  moved  upstream  with  angle  of  attack  and 
that  there  was  attached  flow  over  only  0.17c  for  the  17  degree  angle.  The 
much  higher  values  of  cf  of  the  lower-surface  boundary  layer  are  evident. 


)uas1-Stead<i 


Having  demonstrated  that  the  steady  flow  over  the  Sikorsky  airfoil  can  be 
predicted  accurately  for  all  angles  of  attack  Including  post-stall, 
calculations  were  performed  with  the  quasi-steady  approach.  At  first,  and  to 
check  the  validity  of  the  quasi-steady  model,  two  types  of  slow  motion 
corresponding  to  a  harmonic  oscillation  according  to 


<*(t)  =  5°  1 0°s  1  nu>t  (7) 

-5 

at  a  reduced  frequency  u  of  10  and  a  ramp  type  motion  from  0°  to  16°  at 

-5 

constant  pitch  A(r  da/dt  c/V  )  of  10  were  considered  for  the  same 

00  6 

airfoil  at  a  Reynolds  number  of  2  x  10  .  Under  these  conditions  the  maximum 
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angle  of  attack  exceeded  the  static  stall  angle  In  each  case  and,  since  the 
Influence  of  the  wake  Is  known  to  be  Important,  particular  attention  was 
directed  to  Its  correct  representation.  Subsequently,  calculations  were 
performed  for  higher  frequencies  to  determine  the  unsteady  effects. 

The  variations  of  lift  coefficient,  pressure  coefficients  and  boundary-layer 
parameters  were  examined  for  angles  of  attack  up  to  15  degrees  and  proved  to 
be  Identical  to  those  of  Figs.  2,  3  and  4.  Thus,  the  quasi-  steady-flow 
calculation  method  performed  correctly  and  without  numerical  difficulties. 
Figure  5  shows  the  location  of  the  stagnation  point  and  of  the  displacement 
thickness  at  the  trailing  edge  as  functions  of  time  and,  since  time  Is 
associated  with  angle  of  attack,  the  same  results  could  be  determined  from 
Figs .  2  to  5 . 

The  quasi-steady  calculation  method  was  also  used  to  determine  the  flow  over 
the  airfoil  when  subjected  to  a  ramp-type  motion  from  zero  to  16  degrees  and, 
again,  the  results  were  Identical  to  those  of  Figs.  2  to  4.  Figure  6  Is 
presented  to  show  the  Influence  of  the  wake  on  the  calculation  of  lift 
coefficient  and  confirms  that  It  must  be  Included  for  angles  of  attack  above 
around  8  degrees. 

The  quantification  of  the  wake  effect  In  the  slow-motlon  cases  provides  a 
basis  for  evaluating  the  effect  In  more  severe  unsteady  flows.  When  the 
airfoil  performs  slow  movements,  the  vortlcltles  shed  to  the  wake  are  weak  and 
the  wake,  due  to  viscosity,  behaves  like  a  distribution  of  sinks  along  the 
dividing  line.  Increasing  the  amplitude  and  frequency  of  the  unsteady  motion 
Increases  the  strength  of  the  vortlclty  shed  to  the  wake  and  the  wake  of  the 
flow  over  an  airfoil  can,  therefore,  be  described  as  Including  sources  of 
vortlclty  and  sinks  which  mix  and  Interact  with  each  other.  In  the  present 
calculations,  the  Incorporation  of  the  wake  Involved  the  assumption  that 
vortlclty  and  sinks  are  distinct  and  do  not  Influence  each  other.  This  means 
that  the  vortlclty  shed  via  the  potential  flow  Is  transported  downstream  by 
the  local  velocities  without  the  Influence  of  sinks,  and  the  sinks  computed 
from  the  viscous  wake  are  distributed  on  the  Instantaneous  wake  dividing  line 
whose  location  Is  evaluated  without  the  Influence  of  vortlclty.  This  model 
satisfies  the  steady  flow  condition  as  the  unsteady  motion  approaches  zero  and 
should  capture  most  of  the  characterl sties  of  the  wake  as  the  unsteady  motion 
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of  the  airfoil  Increases.  Some  Instantaneous  wake  dividing  lines  for  the 
Sikorsky  airfoil  executing  a  harmonic  motion,  given  by  Eq.  (7),  at  u  = 

0.001,  0.01,  0.1  and  0.5  are  shown  In  Fig.  7  and  the  wake  locations  for 
frequencies  below  0.1  are  very  close  to  those  for  the  steady-wake.  With  u  = 
0.5,  they  are  quite  different  In  the  far  wake  region,  but  not  In  the  near  wake 
where  the  wake- viscosity  effects  are  most  Important. 

Figure  8  presents  variations  of  lift  coefficient  with  angle  of  attack  for  the 
same  airfoil  subject  to  the  same  ramp  motion  and  same  Reynolds  number  but 
operating  at  higher  pitch  rates  corresponding  to  A  =  0.01  and  0.02.  As 
expected.  Increasing  the  value  of  A  decreases  the  role  of  viscous  effects  so 
that, at  small  angles  of  attack,  the  viscous  lift  coefficient  Is  slightly 
higher  than  Its  Invlscld  counterpart.  The  difference,  however,  diminishes  at 
a  around  9°,  after  which  the  viscous  lift  coefficients  become  less  than  the 
invlscld  values  and  near  and  post-stall  exhibit  the  behavior  of  lift 
coefficients  for  steady  flows.  Figure  8b  shows  that  doubling  the  pitch  rate, 

A  =  0.02,  enhances  the  "over-react"  behavior  of  viscous  lift  coefficients  and 
It  remains  to  determine  If  this  result  Is  real  or  a  consequence  of  the  models 
used  In  our  study.  Comparison  of  Figs.  2,  8a  and  8b  shows  that  an  Increase  In 
pitch  rate  leads  to  an  Increase  In  maximum  lift,  probably  associated  with  a 
delay  In  stall  to  higher  angle  of  attack. 
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4.0  CONCLUSIONS 


The  results  presented  In  this  paper  are  entirely  computational  and  need  to  be 
appraised  In  terms  of  physical  results.  They  show,  however,  that  the 
quasi-steady  computational  method  works  satisfactorily  for  oscillating  and 
ramp-type  motions  of  an  airfoil.  The  Importance  of  the  effects  of  the  viscous 
layer  and  of  the  wake  are  clearly  demonstrated.  In  light  of  the  recent  work 
for  post-stall  steady  flows,  It  would  be  Interesting  and  useful  to  apply  the 
present  method  to  unsteady  flows  over  airfoils  at  angles  of  attack  much  higher 
than  those  considered  here  and  to  Investigate  the  behavior  of  the  solutions  on 
the  airfoil  and  In  the  wake  In  relation  to  measurements.  This  will  Irevltably 
lead  to  the  need  to  replace  the  quas 1 -boundary- layer  approach  by  an 
Interactive  scheme  In  which  unsteady  boundary-layer  equations  are  solved. 
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Fig.  2. 


Variation  of  lift  coefficient  with  angle  of  attack  for  steady  flow 
over  the  Sikorsky  airfoil  at  Rc  =  2  x  106,  —  panel  and  — 
Interactive  methods,  o  data  of  Ref.  18. 


Cakulaed  results  with  the  quasi-steady  model  for  a  harmonic  mot 
according  to  Eq.  (7)  at  u  ,  10-5,  Rc  .  2  x  ,06.  (a)  Locu°  „f  m0t 

5J* tsta9nat,°"  ',01nt’  (H)  variation  of  tralllng-edge 
displacement  thickness  on  the  upper  and  lower  surfaces  with  phase 


Fig.  6.  Effect  of  viscous  wake  on  lift-coefficient  distribution. 

Quasi-steady,  ramp-type  motion  with  a  varying  from  0°  to  16° 
with  constant  pitch  rate,  A  =  10~^,  Rc  =  2  x  10&. 


and 
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a  :  Panel  method 

b  :  Interactive  method  with  viscous  wake 


c  :  Interactive  method  without  viscous  wake 


Fig.  8.  Continued 
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